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We study transverse electron momentum distribution (TEMD) in strong field atomic ionization driven by laser 
pulses with varying ellipticity. We show, both experimentally and theoretically, that the TEMD in the tunneling 
and over the barrier ionization regimes evolves in a qualitatively different way when the ellipticity parameter 
describing polarization state of the driving laser pulse increases. 


PACS numbers: 32.80.Rm 32.80.Fb 42.50.Hz 

A highly non-linear interaction of ultra-short light pulses 
with matter enabled studying electron dynamics on the atomic 
time scale and facilitated emergence of the attosecond sci¬ 
ence [1]. In addition, strong field atomic ionization proved 
itself a potent tool to interrogate atomic and molecular orbital 
structure via high order harmonic radiation [2], tunneling and 
diffraction [3] or tunneling and momentum imaging [4]. This 
utility of strong field atomic ionization is based on the electric 
field of the laser pulse bending the Coulomb barrier and letting 
a bound electron to tunnel out from an atom or a molecule. 

Within the premise of the Keldysh theory [5], this tunnel¬ 
ing regime of strong field ionization corresponds to a small 
value of the adiabaticity parameter (i)^J2Ip/E < 1 de¬ 
fined via the frequency CO and the strength E of the laser field 
and the ionization potential Ip of the target atom (the atomic 
units are used in the paper unless otherwise specified). A 
finer distinction arises when one realizes that the Keldysh the¬ 
ory in its original form is not applicable for very high field 
strengths exceeding the over the barrier (OBI) limit. The OBI 
regime was first observed in [6] (see also [7] for a detailed 
review). The Keldysh theory in its original form fails in the 
OBI regime because there is a classical escape trajectory for 
an electron. One cannot, therefore, rely on the saddle point 
method that Keldysh employed in his original work. The 
so-called Keldysh-Faisal-Reiss (KFR) theory [8, 9] must be 
used instead to describe the OBI regime. We also note that 
the Keldysh theory [5] and its subsequent developments and 
generalizations [10-13] describe the quasistatic limit of small 
laser pulse frequencies. In strict terms, the Keldysh approach 
provides a leading-order term in the asymptotic expansion of 
the ionization rate, a systematic way to obtain higher order 
terms is described in [14, 15]. 

Despite the fact that underlying physics is very different in 
the two regimes (a classically forbidden trajectory for tunnel¬ 
ing and a classically allowed trajectory for OBI), the energy 
spectra and electron angular distributions as given by these 
two theories are not dissimilar. In this Fetter, we demonstrate, 
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that the transverse electron momentum distribution (TFMD) 
is a measurable quantity that is qualitatively different in the 
tunneling and the OBI regimes. This distribution (also known 
in the literature as the lateral electron momentum distribution 

[16] ) gives the probability to detect a photoelectron with a 
given value of the momentum component py perpendicular 
to the polarization plane of the laser radiation. In the tun¬ 
neling regime, TFMD exhibits a cusp-like structure due to 
the Coulomb focusing effect at = 0 for linear polarization 

[17] , and a Gaussian-like structure predicted by the Keldysh 
theory for circular polarization [4]. We studied this transition 
from the cusp-like to the Gaussian structures in detail in the 
tunneling regime [18], and interpreted this transition as a grad¬ 
ual diminishing of the role of the Coulomb effects with grow¬ 
ing ellipticity of the laser pulse. Further study of the role of the 
Coulomb focusing effects was reported in [19]. We shall see 
below that the situation is quite different in the OBI regime, 
where the TFMD always has a cusp regardless of the value 
of the ellipticity parameter. As a result of this qualitatively 
different behavior of the TFMD, one can clearly distinguish 
the tunneling and OBI regimes. This is an important result 
since the TEMD conveys information about the fine details of 
the strong field ionization process [20, 21]. One such detail 
is the electron velocity distributions at the moment of time 
when ionization occurs, which is often used in various models 
of strong field ionization. The omnipresence of the cusp in the 
OBI regime also makes it unsuitable for momentum imaging 
proposed in [4]. 

As the case study, in the present work we select two 
markedly different atomic species: the argon atom in the *5o 
ground state and the neon atom in the metastable state 
with the ionization potentials of 15.76 eV and 5.07 eV, re¬ 
spectively. An estimate for the critical field corresponding to 
the onset of OBI can be found from the equation Fobi — /2 

which follows from considering the hydrogen atom placed in 
an external field using the parabolic coordinates [22]. This 
rough estimate, which does not account for the above-barrier 
reflection [7] and the Stark shift, places the OBI onset of Ar 
well into the 10*^ W/cm^ intensity range while for Ne* this 
onset starts already in the 10^^ W/cm^ range. This com¬ 
fortable two orders of magnitude difference allows to drive 
these targets to the tunneling and OBI regimes with compa- 
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rable laser intensities in the same experimental set up ( Ar@ 
4.8 X 10^"^ W/cm^ and metastable Ne* @ 2 x lO'"^ W/cm^, both 
corresponding to a similar adiabaticity parameter y ~ 0.7 at 
800 nm). To our knowledge the only known momentum imag¬ 
ing experiment in similar OBI regime was reported on Li [23]. 

A schematic representation of the experiment is shown in 
Figure 1. The ultrafast light pulses are produced by a com¬ 
mercially available chirped-pulse amplification laser system 
(Femtolasers, Femtopower Compact Pro CE Phase). The 
light pulses are generated, stretched, amplified and then com¬ 
pressed in the system. The pulse repetition rate is IkHz with 
a pulse duration of 6 fs, pulse energies of approximately 
450 mJ. The pulse train is focused down to a spot size of 
7 pm radius (FWHM) at the interaction region of the Reac¬ 
tion Microscope (REMI). This is the electron detection de¬ 


vice, where the laser pulse ionizes atoms from a target atom 
beam, then uses position dependent delay-line time of flight 
detectors to determine the momentum vectors of the ionized 
electrons. More information on the experimental setup can 
be found in [24]. The electron momentum is measured as a 
function of the ellipticity of the ionizing beam, which is var¬ 
ied using a quarter waveplate. The Ar beam is provided by a 
cold gas jet source. Metastable ^7*2 neon atoms are produced 
by a gas discharge source, which uses a DC discharge across a 
supersonic gas expansion region to excite approximately 1 % 
of neon atoms in a gas jet into the correct state. The flux 
of metastable neon atoms is improved by optical collimation 
techniques that take advantage of the 640 nm closed optical 
transition to the state. Eurther details of this gas source 
can be found in [25, 26]. 
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FIG. 1: Schematic representation of the experiment 


Our theoretical results are obtained by solving the time- 
dependent Schrodinger equation (TDSE); 

= (4,om + Ant(f)) • (1) 

To describe the field-free Ar and metastable Ne* atoms, we 
used effective one-electron potentials [27]. The interaction of 
the atom with the laser pulse is described in the velocity form 
of the interaction operator: 

= A{t) ■ p , (2) 

where A{t) is the vector potential of the laser pulse. The 
laser pulse is elliptically polarized and propagates along the 
z-direction which is assumed to be the quantization axis: 

E Es 

Ex = n -- /(f) COS (Of , Ev = n -y /(0sin(Of , (3) 

VI+e^ ' Vl+e^ 

where e is the ellipticity parameter. The function /(f) = 
sin^(7tf/7j), with T\ being the total pulse duration, is used 


to represent the pulse envelope. Eor the Ar atom, the field 
strength was E =0.1171 a.u. corresponding to the experi¬ 
mental peak intensity of 4.8 x lO^^W/cm^. Eor the metastable 
Ne* atom, E — 0.0756 a.u. with the peak intensity of 2 x 10^^ 
W/cm^. The carrier wavelength A, = 800 nm and the EWHM 
of 6 fs were the same for Ar and Ne*. 

To solve the TDSE we employ the strategy used in the pre¬ 
vious works [18, 28, 29]. The solution of the TDSE is repre¬ 
sented as a partial waves series: 

Enil x I 

1=0 fj=-l 

The radial part of the TDSE is discretized on the grid with the 
stepsize 5r = 0.1 a.u. in a box of the size Rmax = 400 a.u. 
The maximum orbital momentum in Eq. (4) was restricted 
to Emax = 60. Convergence with respect to variation of 5r, 
^max and Lynax was Carefully monitored. The matrix iteration 
method [30] was used to propagate TDSE in time. Ionization 
amplitudes a{p) were obtained by projecting the solution of 
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FIG. 2: (Color online) TEMD of Ar (multiplied by 10^) for ellipticity 
parameters e = 0, 0.42, and 1 (from top to bottom). TDSE calculation 
is shown by a (red) solid curve (shaded for a clearer appearance), 
experimental data are plotted with error bars. The peak intensity of 
the laser pulses is 4.8 x lO*'^ W/cm^. 


FIG. 3: (Color online) TEMD of metastable Ne* for ellipticity pa¬ 
rameters e = 0, 0.42, and 1 (top to bottom). TDSE calculation is 
shown by a (blue) solid curve (shaded for a clearer appearance), ex¬ 
perimental data are plotted with error bars. The peak intensity of the 
laser pulses is 2 x 10*^ W/cm^. 


the TDSE after the end of the pulse on the set of the ingoing 
scattering states 'Ep (r) of the target atom. The TEMD W (p±) 
was calculated as 

W{p_l) = j \a{p)\^ dpx dpy , p±=Pz (5) 

Experimental and theoretical TEMD results for Ar are 
shown in Eigure 2. A general trend of the calculated TEMD 
with increase of ellipticity is very similar to that reported pre¬ 
viously for the simulations of the hydrogen atom [18]. The 
cusp-like structure is present for linear polarization and it 
gradually evolves into a Gaussian distribution as the elliptic¬ 
ity parameter increases. Agreement between the theory and 
experiment is good for linear polarization but gradually dete¬ 
riorate with an increase of ellipticity. 

Eigure 3 presents the theoretical and experimental TEMD 
results for metastable Ne*. In this target atom, the TEMD 
evolution with the ellipticity parameter is greatly reduced with 
the cusp clearly present even for the circularly polarized pulse. 
Similar to Ar, agreement between the theory and experiment 
progressively worsens from the top to bottom panels. 


To analyze the cusp in more detail, we zoomed in on the 
narrow range of momenta \pi_ \ < 0.25 a.u. and analyzed the 
function V{pi_) = \nW{p\^ in this interval. Eor the TEMD 
VT(px) to have a cusp, V{p±) should have an infinite deriva¬ 
tive of some order and have an expansion near = 0; 

E(px)=B+A|pxr. (6) 

Such expansions, in fact, reproduce V (px — 0) very well as 
was shown in [18]. 

The same functional form (6) was used to fit both the the¬ 
oretical and experimental data for the ground state Ar and the 
metastable Ne* in the whole range of ellipticities by consid¬ 
ering the coefficients A, B, a as fitting parameters. The most 
essential a parameters are shown in Eigure 4 for Ar (top) and 
Ne* (bottom). Both theoretical and experimental values are 
shown with error bars which represent the fitting error. Eor 
the theoretical data, this error does not exceed a fraction of a 
percent and is not visible on the scale of the figure. 

The a parameters shown on the top and bottom panels 
of Eigure 4 demonstrate a qualitatively different behavior as 
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FIG. 4: (Color online) The fitting parameter a in Eq. (6) as a function 
of the ellipticity parameter e for Ar (top) and Ne* (bottom). TDSE 
results are shown with squares (a smooth solid line is to guide the 
eye). The experimental data points are shown with error bars. 


functions of the ellipticity. For the Ar atom, the a parameter 
grows with e reaching the value close to 2 for circular polar¬ 
ization. This implies that TEMD becomes close to a Gaussian 
W{pi_) exp(—with the Gaussian width related to 
the fitting parameter (3—1/ \J—A . The corresponding numer¬ 
ical values of 0.25 ± 0.002 and 0.28 ± 0.02 for the TDSE and 
experiment, respectively, are close to the experimental values 
reported in [4] for comparable field intensities. In the mean¬ 
time, the a parameters for the metastable Ne* atom remain 
essentially flat, indicating that a cusp-like behavior is present 
for all e in the range from linear to circular polarization. In this 
case, extraction of the Gaussian width parameter is not possi¬ 
ble even for the circular polarization. Several TDSE calcula¬ 
tions performed for different field intensities did not show any 
considerable variation of the cusp width. However, the Gaus¬ 
sian width varies with the field strength as o This may 

explain, at least partially, deviation between the measured and 
calculated TEMD with circular polarization due to the varia¬ 
tion of the field strength across the laser-atom interaction re¬ 
gion while the calculation was performed at a single nominal 
field intensity. 

While the Ar case shows the behavior qualitatively simi¬ 
lar to that found previously for hydrogen [18], the metastable 
Ne* presents a different trend, with the cusp never disappear¬ 
ing completely. In this case, a simplified description based on 
the Keldysh theory is never correct even qualitatively. This 
qualitative difference can be explained by the different ioniza¬ 
tion regimes for Ar and metastable Ne*. 

The TEMD cusp disappearance with increasing e can be 
related to a dramatic change of the angular momentum com¬ 
position of the ionized electron wave function [18]. This com¬ 
position is characterized by the distribution of the norm Ni of 
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FIG. 5: (Color online) Angular momentum distribution A; for Ar@ 
4.8 X lO^'* W/cm^ (top) and metastable Ne*@ 2 x 10^'^ W/cm^, bot¬ 
tom). Laser field is circularly polarized. 


the wave function obtained if only the terms with spherical 
harmonics of rank I are retained in expansion (4). For a tun¬ 
neling process this distribution is shifted towards larger I with 
increasing pulse ellipticity parameter. Indeed, tunneling can 
be viewed as a non-resonant absorption of a large number of 
photons. Absorption of a photon from the circularly polarized 
wave increases the magnetic quantum number by one unit. 
This leads to a prevalence of high angular momenta in the par¬ 
tial wave expansion (4). High angular momenta create large 
centrifugal barrier preventing recolliding electron trajectories, 
thereby suppressing the Coulomb focusing effects. The cusp, 
therefore, vanishes for polarization close to circular, as in the 
case of Ar reported here, or hydrogen [18]. The situation with 
the metastable Ne* is completely different. OBI dominates 
in this case, and since OBI is essentially a distortion of the 
atomic potential to the degree, that there is effectively a zero 
barrier to the continuum, the atom does not have to absorb 
many photons to become ionized. The distribution A/, there¬ 
fore, is peaked at lower values of the angular momenta. That 
this is indeed the case can be seen in Figure 5, where the dis¬ 
tributions A/ are shown for Ar for the intensity of 4.8 x 10*"*^ 
W/cm^ (tunneling) and Ne* for the field intensity of 2 x 10^^ 
W/cm^ (the OBI regime). Smaller angular momenta enhance 
the area near the origin where the Coulomb focusing effect 
is strongest. Larger angular momentum components are re¬ 
pelled from the origin due to the centrifugal barrier. Hence 
in the former case the cusp is always present whereas in the 
latter case it gradually vanishes. This corresponds to a clas¬ 
sical trajectory starting from the origin whereas a tunneling 
trajectory starts at the point of exit from the tunnel. In the 
OBI regime the electron’s classical trajectory starts at the ion 
core regardless of the polarization of the laser pulse, which 
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may be enough for the efficient Coulomb focusing even if the 
trajectory never returns to the core. 

To summarize, we described an effect of bending the 
Coulomb barrier of the atom on the transverse electron mo¬ 
mentum distribution (TEMD) in strong field ionization in the 
tunneling regime. This fundamental effect, which should be 
present in any atomic or molecular target, is measured experi¬ 
mentally and modeled theoretically in two markedly different 
species: the ground state Ar and metastable Ne*. The effect 
is substantial, it has never been described or observed before 
and it enables a clear distinction between the tunneling and 
OBI regimes in the experiment. Also, it has to be taken into 
account when using TEMD data to interrogate electronic or¬ 
bitals of the target. 


Einally, we note that the cusp disappearance in the case of 
circular polarization may seem to follow from a classical con¬ 
sideration. Indeed, in the circularly polarized field, the two or¬ 
thogonal field components drive the photoelectron away from 
the ionized core thus reducing the Coulomb focusing effect. 
This classical consideration, however, fails to explain quali¬ 
tatively different TEMD behavior in the tunneling and OBI 
regimes observed in the present study. 
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